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Low-temperature transport spectroscopy measurements on a suspended few-hole carbon nanotube 
quantum dot are presented, showing a gate-dependent harmonic excitation spectrum which, strik- 
ingly, occurs in the Coulomb blockade regime. The quantized excitation energy corresponds to the 
scale expected for longitudinal vibrations of the nanotube. The electronic transport processes are 
identified as cotunnel-assisted sequential tunneling, resulting from non-equilibrium occupation of 
the mechanical mode. They appear only above a high-bias threshold at the scale of electronic nan- 
otube excitations. We discuss models for the pumping process that explain the enhancement of the 
non-equilibrium occupation and show that it is connected to a subtle interplay between electronic 
and vibrational degrees of freedom. 

PACS numbers: 63.22.Gh, 73.63.Fg, 73.23.Hk 



The coupling of vibrational modes and electronic 
transport in nanoscale systems and in particular quan- 
tum dots is currently the focus of many theoretical 
[H [2] and experimental [3 H E O H H [9] research 
efforts. In this respect, single wall carbon nanotubes 
(SW-CNTs) provide a unique mesoscopic system where 
both bulk beam mechanics [3l [4] and quantization of 
phonon modes [5 , 6, 9 have already been demonstrated. 
Low-temperature Coulomb blockade (CB) spectroscopy 
on quantum dots formed within the nanotube [lOl [TT] 
has led to a well-developed understanding of the elec- 
tronic structure of SW-CNTs [12| [131 [M] . suspended 
SW-CNT quantum dots, transport spectroscopy has also 
revealed the Franck-Condon effect [T, where the quan- 
tized vibrations of the nanotube become visible in single 
electron tunneling (SET) at finite bias. A large electron 
phonon coupling [6 , 9 , and first indications of vibrational 
phenomena in cotunneling have been observed [9^ . 

In this Letter we present low-temperature transport 
measurements on a suspended CNT quantum dot sys- 
tem. We observe signatures of non-equilibrium popula- 
tion of the quantized mechanical oscillations in the trans- 
port spectrum, revealed by cotunnel-assisted sequential 
electron tunneling (CO-SET). The non-equilibrium oc- 
cupation is enhanced ( "pumped" ) by higher-order tunnel 
processes. Detailed models are discussed which explain 
the observations. 

The basic device geometry is sketched in Fig. [TJa). A 
highly P++ doped silicon wafer, also functioning as back 
gate, with 500 nm thermally grown oxide on top provides 
the starting point of the sample preparation. After litho- 
graphic fabrication of marker structures and localized de- 
position of growth catalyst [16 , carbon nanotubes are 
grown in situ by chemical vapor deposition and located 
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FIG. 1: (Color online) (a) Side view drawing of the sample 
geometry, (b) Overview measurement of the differential con- 
ductance d//dVsD of the nanotube quantum dot system in 
the few- hole region (1 < A^h < 2, with A^h = and A^h = 3 
visible at the edges of the plot) as function of back gate volt- 
age Vg and bias voltage Vsd (logarithmic color scale, negative 
differential conductance is plotted white). In the CB regions, 
the number of trapped valence band holes A^h is indicated. 



using atomic force microscopy. Contact electrodes con- 
sisting of 5 nm chromium and 50 nm gold are deposited 
on top of the nanotubes. In a last step, the nanotube 
devices are suspended using wet etching in buffered hy- 
drofluoric acid. 

Measurements were performed in a dilution refrigera- 
tor with a base temperature Tmc ^ 20 mK and an elec- 
tron temperature Tq\ 100 mK. Fig. [TJb) shows the dif- 
ferential conductance dJ/dFsD as function of back gate 
(substrate) voltage Vg and bias voltage Fsd of one partic- 
ular suspended nanotube quantum dot. This sample has 
a lithographically designed length of / = 250 nm, and was 
- from transport measurements as in Fig. [ijb) - identi- 
fied to be a semiconducting nanotube with a band gap of 
Eg = 200 meV. Band gap energy and device length lead 
to a predicted energy scale of ^ 6 meV of orbital elec- 
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tronic excitations [14]. In the measurement of Fig. [TJb), 
the basic structure of diamond-shaped CB regions with 
a fixed trapped charge, as expected for a single quan- 
tum dot, is clearly visible. Since the band gap position is 
known, we can identify the charge states as A/'h = 1 and 
A/h = 2, respectively. 

Figure [2|a) displays a detail measurement of the sin- 
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FIG. 2: (Color online) (a) Differential conductance d//dVsD 
in the region where the nanotube is charged with 1 < A^h < 2 
holes, as function of gate voltage Vg and bias voltage Vsd (lin- 
ear color scale), (b) Detail zoom of part of the SET region 
with 1 < A^h < 2, as marked in (a), using a different color 
scale. In both (a) and (b), arrows point out line features cor- 
responding to excited states (see text) [17]. (c) Excitation 
energies corresponding to the lines of enhanced d//dVsD in 
(a) and (b) with positive slope; the x-axis is the line number. 
The solid line is a linear fit, resulting in an average energy dif- 
ference of = 0.425 it 0.004 meV per line, (d) Detail zoom 
of the Coulomb blockade region with A^h = 2, as marked in 
(a) (logarithmic color scale) [17]. A white dashed line sketches 
the edge of the SET region. For better contrast, different color 
scale ranges are chosen in parts of the plot, (e) Relative ex- 
citation energies for the line features in (d), using the same 
parameters for position to energy conversion. A linear fit re- 
sults in A^; = 0.810 ± 0.025 meV. (f) Schematic detaiUng 
cot unnel- assisted sequential tunneling (CO-SET) (see text). 

gle electron tunneling region with 1 < A^h < 2 at low 
negative bias, and the adjacent A/'h = 2 CB region. A 
rich spectrum of equidistant excitation lines with posi- 
tive slope, corresponding to excitations of the A^h = 1 
system (see e.g. Ref. [14 ), is found in SET [17 . This 
is also detailed in the enlarged plot of Fig. |2|b) (see ar- 
rows). Figure |2|c) shows the corresponding excitation 
energies as function of line number ^21^ . We assign these 



excitations to a harmonic vibration whose frequency is 
hij = 0.425 meV ± 0.004 meV, in good agreement with 
the bulk mechanics prediction of ficc^vib = 0.44 meV for 
the longitudinal vibration mode of a 250 nm long nan- 
otube segment [6^. In addition, the data of Fig. |2|a-b) 
reveals three faint excitation lines with negative slope 
(i.e. A^h = 2 excitations), marked by three black arrows 
and separated by approximately 0.7 meV. It is difficult 
to confirm a harmonic spectrum because of the faintness 
of the three lines. 

Figure |2|d) enlarges the region outlined in Fig.[2|a) by 
a black dashed rectangle, plotting the differential con- 
ductance in logarithmic color scale. Here, CB stabi- 
lizes a total charge of A^h = 2 holes on the nanotube, 
suppressing SET. Surprisingly, a pattern of gate depen- 
dent excitation lines in the CB region, parallel to the 
edge of the SET region, emerges [IT. Their relative 
excitation energies are plotted in Fig. [2|e). A regular 
spacing corresponding to a harmonic oscillator energy of 
huj = 0.810 meV±0. 025 meV is visible, close to the above 
mentioned energy scale 0.7 meV of the A/'h = 2 excitations 
in SET. 

The most straightforward explanation for these lines 
is that they correspond to cotunnel-assisted sequential 
tunneling (CO-SET) processes [18 . These are multi-step 
processes, as sketched in Fig. [2|f). Inelastic cotunneling 
leaves a quantum dot in an excited state (see second panel 
in the figure), which is possible at any voltage above the 
energy of this state. Subsequently, sequential tunneling 
to the drain electrode can follow for voltages close to the 
Coulomb diamond edges, even if the Coulomb-blockaded 
ground state does not allow it. For this to happen, the 
rate F for tunneling out has to be comparable to or larger 
than the energy relaxation rate 7 to the ground state. 

In general, CO-SET processes can involve intrinsic 
(e.g. orbital) excitations of quantum dots as well as vi- 
brational states. Here, as opposed to the measurements 
of Ref. [18 , the highly regular excitation spectrum ob- 
served in Fig.[2|d) strongly indicates a vibrational origin. 
The vibrational CO- SET lines are related to the two- 
hole charge state vibrational excitations seen in SET with 
negative slope. As the vibrational energy exceeds ther- 
mal broadening {hw ^ k^T)^ the observation of multiple 
excitations in CB strongly suggests storage and subse- 
quent release of energy in the vibrational mode, involving 
phonon absorption sidebands [2] [S] [19] . This stands in 
contrast to the usual assumption that vibrational relax- 
ation is fast and that the mechanical system is predom- 
inantly found in its ground state. The required energy 
is provided by inelastic cotunneling, exciting the vibra- 
tional mode and causing its non-equilibrium occupation. 

A particularly interesting feature of the measurement 
of Fig. |2|a-e) is that the CO-SET lines can only be 
observed beyond a weakly gate-dependent threshold at 
— — 4mV (marked X in the plots), where the dif- 
ferential conductance is enhanced; i.e. they only occur 
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at energies higher than 4 meV. The hne feature marked 
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FIG. 3: (Color online) Single quantum dot model reproduc- 
ing the observed transport features: (a) Process and energy 
level scheme analogous to Fig. |2jf), leftmost panel, show- 
ing the lead Fermi levels, the accessible states within the 
quantum dot, and tunneling and relaxation rates relevant for 
pumping the vibration mode into non-equilibrium [17 . For 
a detailed description of the panels and the rates F, F^, and 
7^ see the text, (b) Calculated differential conductance as 
function of gate voltage Vg and bias voltage Vsd, for the fol- 
lowing parameters: vibrational level-spacing hou = 810/ieV, 
kT = 8.6/ieV = lOhF' = lO^hF = lO'^hy = 10^/17. The tun- 
nel couplings were chosen smaller than in the experiment to 
ensure a well-behaved perturbation expansion. Vg is scaled 
with the gate conversion factor a = Cg/C (see e.g. [151 fTSl ). 
Arrows indicate the enhancement of CO-SET at high bias. 

with X in the measurements represents the onset of an 
inelastic cotunneling current, corresponding to an elec- 
tronic excitation of the quantum dot. To understand 
the resulting interplay between mechanical and electronic 
excitations we have performed model calculations, illus- 
trated in Fig. [3) The transport calculations account for 
strong Coulomb effects as well as tunneling processes up 
to 4th order in the tunneling Hamiltonian, responsible 
for cotunneling, line broadening and level renormaliza- 
tion which are important close to SET resonance; see 
Ref. [19| for details. 

The two-hole ground state and its vibrational excita- 
tions are the origin of CO-SET lines as observed in the 
experiment. To demonstrate this idea, it suffices to con- 
sider a model with several equidistant states A-E (see 
Fig. [sj^a), panel (1), and also [IZ]), which are coupled 
with the same rate F to the one- hole ground state. In 
addition, an excited two-hole state X supporting a cur- 
rent with rate F' is introduced. Since the edge of the 
CB diamond becomes very broad in the experiment be- 
yond this excitation, an estimate of the ratio of F and 
F' based on the SET current is difficult. However, the 
strong signature of the inelastic cotunneling line associ- 
ated with state X indicates that F' F. The result of 
a calculation based on this model is shown in Fig. [sj^b). 
The important characteristic of the experiment is repro- 
duced: CO-SET lines are strongly enhanced beyond the 



threshold for inelastic cotunneling connected to state X 
(see arrows in Fig. [sjb)). 

The appearance of a set of harmonic CO- SET lines 
at high bias can be explained by a subtle interplay of 
electronic and vibrational excitations, with two coexist- 
ing paths. Figure [sja) illustrates these two paths; panels 
(1) and (2) for one and panel (1) and (2') for the other. 
For energies larger than the one of the inelastic cotunnel- 
ing step at VsB — —4 mV a process involving two stages 
of inelastic cotunneling becomes feasible: the first effec- 
tively excites the dot from state A to state X (panel (1)), 
the second results in a transition from state X into one 
of the states B-E (panel (2)). The rate for each such 
cotunnel process is oc FF^ Compared to the low bias sit- 
uation, this significantly enhances the population of the 
vibrational excited states B-E, and thereby the visibil- 
ity of the vibrational CO-SET lines, since without access 
to state X, the limiting rates for CO-SET are propor- 
tional to F^ « FF'. Thus, while CO-SET lines can in 
principle also be present at small bias (|1^sd| < 4mV), 
beyond the inelastic contunneling step they appear more 
pronounced: The strongly coupled excited state pumps 
the vibrational mode out of equilibrium, enhancing its 
population. 

The second path that increases the occupation of the 
vibrationally excited states B-E involves direct relaxation 
of state X (panel (2') in Fig.jsja)) with a rate 7^ As the 
precise value of 7' is unknown, we have verified that the 
qualitative result that the CO- SET lines are strongly en- 
hanced beyond the inelastic cotunneling step connected 
to state X persists for a large range of this parameter. An 
important further qualitative conclusion from the obser- 
vation of CO-SET features is that relaxation from states 
B-E into the two- hole ground state A (rate 7 in Fig. |2|f)) 
in the experiment can not be faster than tunneling out 
of the quantum dot (7 < F). 



i(a) (b) 




FIG. 4: (Color online) Alternative mechanism for obtain- 
ing a vibrational non-equilibrium occupation, (a) Schematic 
side-view drawing detailing the possible formation of a small 
quantum dot beneath one of the metallic contact leads, (b) 
Transport processes for this 'double quantum dot' case: addi- 
tional cotunneling processes are enabled when the small quan- 
tum dot enters the energy window given by Vsd (see text). 

So far we have assumed that the weakly gate- 
dependent feature X in Fig. [2] corresponds to an elec- 
tronic excited state. An alternative explanation invokes 
an electronic ground state of a second quantum dot form- 
ing within the nanotube, e.g. beneath one of the contact 
electrodes, as drawn in Fig. [4]^ a). Assuming this state has 
a much weaker gate coupling, we arrive at an equivalent 
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mechanism for the non-equihbrium occupation of the vi- 
brational mode. Local distortions and curving edges in 
the measurement of Figs, [l] and [2] may be taken as a hint 
of such a scenario, as has also been reported in Ref. [20] . 

Fig.Qb) illustrates the transport process for this 'dou- 
ble quantum dot' case. The shielded additional small 
quantum dot displays a much weaker gate voltage de- 
pendence of the level energy. The line of enhanced dif- 
ferential conductance marked in Fig. [2] with X now de- 
scribes its alignment at the drain lead with the drain 
Fermi edge. The small quantum dot can then be oc- 
cupied by a sequential tunneling process as sketched in 
Fig. ^h). Inelastic cotunneling through the main dot 
may follow, leading to the occupation of an excited vi- 
brational state. Although the origin of, and tunneling 
processes associated with, state X differ from the first 
scenario, the main effect is the same: state X pumps the 
vibrational mode through consecutive tunneling events. 

An intriguing feature of the measurement that has not 
been discussed so far is that the vibrational frequency in 
the A/'h = 2 charge state, as measured both from SET and 
CO-SET lines, is approximately twice as large as that in 
the A/'h = 1 charge state, seen as SET features only. This 
difference does not affect the pumping mechanism. It is 
possibly related to the fact that the observed quantum 
dot is in the few-carrier regime - in contrast to previ- 
ous work on longitudinal phonon excitations in carbon 
nanotubes |6]. Whereas tension only affects the longi- 
tudinal mode via higher-order effects [6 , the transition 
from e.g. one to two trapped holes involves a distinct 
spatial redistribution of charge along the nanotube [12] . 
This strongly affects the electrostatic force distribution, 
which makes variations of vibration mode energies and 
mode shapes likely. As shown in Ref. |6| , the coupling to 
the lowest symmetric vibrational mode of the nanotube 
may vanish for a charge distribution localized at the cen- 
ter of the nanotube. In general one may thus expect that 
the excitation of different modes depends strongly on the 
charge state, as we observe in the experiment. 

Finally, the experimental observation of vibrational 
CO-SET resonances allows to establish a lower bound- 
ary for the quality factor of the longitudinal mechanical 
mode [2 , independent of the detailed pumping mecha- 
nism. The tunnel current at the edge of the low-bias 
CB region (|Vsd| < 4mV), / 1 nA, provides an ap- 
proximation for the SET tunneling rate F. For CO-SET 
to be visible, the lifetime of vibrational excitations must 
be larger than the corresponding timescale r 0.16 ns. 
With the vibrational energy hco = 810 /ieV, one obtains 
Q > r X u;/27T = 31. Since we observe several sidebands, 
a higher value for Q is more likely than this lower bound- 
ary, depending on the specific relaxation processes. 

In conclusion, we have observed a pattern of equidis- 
tant, gate-dependent vibrational excitations in the 
Coulomb blockade regime of a suspended carbon nan- 
otube quantum dot. The appeareance of the excita- 



tions is explained in the context of co-tunnel assisted 
sequential electron tunneling via phonon absorption pro- 
cesses. Interestingly, the absorption sidebands of the 
quantized phonon mode are only visible above a finite 
bias voltage threshold. Two models for the enhance- 
ment of the non-equilibrium distribution of the mechani- 
cal mode are discussed which demonstrate that the mode 
can be "pumped" either by an electronic excitation of the 
nanotube, or by a ground-state of another small quantum 
dot within the same nanotube. The pumping mechanism 
offers a perspective for electric control of quantized me- 
chanical motion in nanoscale transistors by employing 
the interplay with electronic degrees of freedom. 
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